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Abstract 

We study two-dimensional classically integrable field theory with independent boundary con¬ 
dition on each end, and obtain three possible generating functions for integrals of motion when 
this model is an ultralocal one. Classically integrable boundary condition can be found in solving 
boundary K± equations. In quantum case, we also find that unitarity condition of quantum R- 
matrix is sufficient to construct commutative quantities with boundary, and its reflection equations 
are obtained. 


1 lutroductiou 


Recently, there has been great progress in understanding two-dimensional integrable field theory on a finite 
interval with independent boundary condition on each end 0-i- The motivation is not only from the 
necessary of itself, but also from the studies in boundary-related phenomena in statistical systems near 
criticality and integrable deformations of conformal field theories [Q. 

In order to deal with integrable models with boundary, relying on previous results of Cherednik 0, 
Sklyanin 0, 0 introduced a new generating function which originates from the periodic boundary one. In 
classically integrable models j^, if there has a well-known relation for monodramy matrix 0 as {TfR} = 
[r, T 0 T], and r-matrix satisfies the condition of r{a) = —r(—a), the new generating function defined in 
[x_, a;+] can be expressed as: 

T(a) = tr {K+(a)T{x+, X-, t, a)K-{a)T~^{x+, a;_, t, —a)} , (1) 

where K± are boundary reflection matrices. 

Expanded as a Laurent series in a, all coefficients of T{a) make an infinite number of integrals of motion 
which ensure the completely integrability of the model. From Q, T{a) must be in involution between 
different spectral parameter a and it is independent of time. In other words, K± must satisfy some constraint 
equations, and existence of nontrivial K± solutions means there are nontrival classically integrable boundary 
conditions (CIBC). 

There has no condition of r{a) = —r{—a) in affine Toda field theory (ATFT), so P.Bowcock et al Q 
developed a method of modified Lax pair to deal with such models, in which the new generating function in 
(—oo,a;+] reads 

T{a) = tr {T^(—oo, a:+, t, —a)Kj^{a)T{—(X), a;+, t, a)} , (2) 

in which ”f” denotes conjugation and it has little difference with the original paper0| according to different 
definition in T matrix. We must point out that boundary Lax pair in 0 has been modified from the periodic 
boundary one. 

In this paper, we find it is necessary to add a new parameter to the generating function (0) in order to 
deal with ATFT, and no symmetry condition of r- matrix is needed in fact. Besides this modified form, we 
also construct two other possible generating functions by zero curvature representation. After we extend 
our results to quantum integrable systems, we find unitarity condition of quantum R- matrix is sufficient to 
construct commutative quantities with boundary too. 

The paper is organized as follows. In section 2, three possible generating functions are constructed by 
zero curvature representation. In order to regard constructed functions as generating functions, algebra 
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equations (reflection equations) and evolution equations of K± matrices appear in section 3. Then, we study 
ATFT in section 4 and find the links among these generating functions. In section 5, we extend our results to 
the quantum case, and demonstrate that unitarity condition of quantum R- matrix is sufhcient to construct 
commutative quantities. Then, we compare our commutative quantities with those of paper |l3| and find 
the relation between them. At last, some discussion will be found in section 6. 


2 Construction of generating function 

2.1 Periodic boundary condition 

The zero curvature approach to inverse scattering relies on the existence of a pair of linear partial 
differential equations in d x d matrix. 

dx'^ = U{x,t,a)'^, dt'^ = V{x,t,a)"^, 

here. Lax pair {U,V) are d x d matrices whose elements are functions of complex valued field (f>{x,t) and 
its derivatives, a G C is a spectral parameter. Zero curvature condition appears from compatibility of the 
above equation, it is 

dtU - dxV + [U, P] = 0. (3) 

By zero curvature representation, we define transition matrix 

T{x,y,t, a) = V exp U{x',t,a)dx'^ , x>y, (4) 

where V denotes a path ordering of non-commuting factors. Now, T matrix satisfies 

dxT = U{x,t,a)T, 

dtT = V{x,t,a)T-TV{y,t,a), (5) 

Id = T{x,x,t,a), 


where Id is d x d identity matrix. 

It is well known that trace of monodramy matrix Ti(t, a) = T{L,—L,t,a) is a generating function on 
periodic boundary condition, so is another more explicit form r(a) = hitrTL{t, a). Expanded as a Laurent 
series in a, T{a) make an infinite number of integrals of motion. Conservation condition of these integrals 
can be proved by the second equation of (|^) with periodic boundary condition, and the involution condition 
is proved in Poisson bracket: 

{T{x,y,a)®T{x,y,P)} = [r{a,P),T{x,y,a)'SiT{x,y,0)] ,L>x>y>-L, ( 6 ) 
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in which, T is a d x d matrix, T= T®Id and T= Id®T. r{a, (3) is a d^ x d^ matrix whose elements depend 
on a and (3 only. The Jacobi identity for the bracket holds if and only if r-matrix is a solution of classical 
Yang-Baxter equation. 

In periodic boundary condition, it is obvious that {r(Q:) ® t(/ 3)} = 0. So T{a) constructs a family of 
generating function for integrals of motion. 

2.2 Independent boundary condition 

As soon as periodic boundary condition is broken, T{a) defined before should be not a conservative quantity, 
so that we have to find a new expression of generating function. As discussed in papers |^, if Lagrangian 
density in bulk theory is £/, then, the new Lagrangian density with boundary appears as: 

C = 6{x+ - x)9[x - X-)Cf - 6{x+ - x)V+{4>[x+),d^4>{x+)) - S(x - x-)V-(^(x-), d^^(x-)). (7) 

By means of principle of the least action associated with (|^), we will obtain motion equation in (x-,x+) 
and boundary equations on each end. 


2 


If £/ is expressed as £/ = \d^(j)d^(j) — V{(j)) and V± depend only on (j){x±) (and independent of its 
derivatives), boundary equations on each end will be 


dx4' = Td 4 ,V± , X = x±. 


( 8 ) 


On each end, if (/>(a:) is a smooth function of coordinate x, we can extend the motion equation to 
whole domain [a:_,a;+], so are the boundary Lax pair. In other words, we keep the uniform expression of 
fundamental Poisson bracket in the whole domain, ever if it is on independent boundary condition. 

Similar to the definition of transition matrix, there is another matrix function 


F{x, t, a) 


V 



t>to. 


By zero curvature condition, we construct a quantity which is independent of time: 


(9) 


F ^(a;+, ti, a)T{x+,X-,ti,a)F{x-,ti, a) 


F ^(a;+, t 2 ,a)T{x+,X-,t 2 , a)F{x-,t 2 , a). 


( 10 ) 


It can be proved easily because both sides in the equation are equal to T(a;+, X-,to,a). 

With another equation of argument {—a + S), which can be obtained by the same method, we expect 
generating function with boundary shall be constructed as follows. For example: 

tr { [F{x+,t,—a + 6)F~^{x+,t,a)] T{x+,X-,t,a) 

X [F(x-,t, a)F~^(x-, t, —a + 5)] T“^(a;+, a;_, t, —a + (5)} . 

By its product method, it is obvious that this quantity is a conservative quantity. If we regard it as a 
generating function for integrals of motion, there are the main problem: 1. Does such a constructed quantity 
satisfy involution condition ? 2. It had better be independent of ”to” which comes from F matrix. Base 
on these problems, we introduce K± matrices instead of F terms and impose involution and conservation 
conditions on the new form. The new generating function is 

a. 

T{a) = tr {iF+(a:+, t, a)T{x+,X-,t, a)K-{x-,t, a)r“^(x+, a;_, t, —a + (5)} . (11) 

Moreover, we use ’’f” (conjugation) or (transposition) instead of ” —1” (inverse) in order that iF+ 
and K_ depend only on the variables of the boundary x+ and a:_, respectively. The results are: 

b. 

r(a) = tr {iF+(a;+, t, a)T{x+, X-,t, a)K_{x-,t, a)T*{x+, a;_, t, —a + 5)} , (12) 

c. 

r(a) = tr {iF+(x+, t, a)T{x+, X-,t, a)K-{x-,t, a)T^{x+,x-,t, —a + ^)} . (13) 

In quantities and (0), K± matrices and 6 are similar in symbols only. Each of these quantities will 

be a generating function of integrable systems with boundary, if both involution and conservation conditions 
are satisfied. 


3 t(q;) as a generating function 


It is well known that generating function for integrals of motion must be in involution (between each other) 
and independent of time |^. So if we regard quantity & as a generating function, some constraint conditions 
must be imposed on it. Now we study quantity & in this section. 

3.1 Involution condition 

Taking the notations similar to paper [^, we define: 

T+{x,a) = T~^{x+,x,t,—a + S)K+{x+,t,a)T{x+,x,t,a), 

T_{x,a) = T{x,X-,t,a)K_{x-,t,a)T~^{x,X-,t,—a + S), 

T{x,a) = T_(x,a)T^{x,a). 
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Comparing it with quantity (0, we find trT{x,a) is just equal to T(a). Now, we impose some constraint 
conditions on K± matrices: 

[K±{x±,t,a) f T{x+,x-,t,l3)) = 0, 

{X±(a)®if±(/3)} = 0 , {if±(a)®i^^(/3)} =0. (14) 

It means 

{T+{x,t) f T_(x,t)} = 0. 

If K± matrices are independent of field variances, condition is satisfied naturally. But we must point 
out that the K± matrices, in general, can depend on the field variables. Poisson bracket on 7± is: 

Proposition 1. If K+ matrix satisfies 

0 = -r{-a + 5,-^ + 6)K+{t,a)K+{t,'f) + K+{t,a)r{a,-^ + S)K+{t,^) 

+K+it, 7 )r(-a + 5, j)K+{t, a) - K+{t, a)K+it, j)r{a, 7 ), (15) 

then, T+ algebra should obey the following relation: 

|t + {x,a),T+ (a:, 7 )| 

12 1 2 
= -r{-a + (5, -7 + S)T +(a:, a)T +ix,j) + T +{x,a)r{a, -7 + S)T +{x,^) 

+T +(a:, 7 )r(-a + ( 5 , 7 )r +(a:,Q:) -T+{x,a)T +(x, 7 )r(a, 7 ). (16) 

It should be emphasized that K+ is a subalgebra of T+ algebra according to the definition of 7i|_. Propo¬ 
sition 1 can be proved by calculating Poisson bracket on 7i|_ directly. There is another algebra of 7T similar 
to 7i)_: 


Proposition 2. If K- matrix satisfies: 

0 = r{a,j)K-{t,a)K-{t,^)-K-{t,a)r{-a + S,j)K-{t,^) (17) 

- K- {t, 7 )r(a, -7 -f S) K- {t, a)-l- K- (t, a) K- {t, j)r{-a -I- <5, -7 -f <5) , 
then, it leads to the relation of 7T algebra being: 

|r- {x,a),T- (a:, 7 )| 

12 1 2 

= r{a,^) T- {x,a) T- (x,^)-T- {x,a)r{-a + S^'f) T-(x,^) (18) 
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- T- {x, 7 )r(Q!, -7 -I- S) T- (x, a)+ T- (x, a) T- (x, 7 )r(-a -f <5, -7 -f <5). 


The proof is similar to proposition 1. Used proposition 1 and 2, Poisson bracket on T(x,a;) can be 
calculated as follows: 

r (x,a),r (x, 7 ) 

= |r- (x,Q!) T+ (x,a),T- (x,7) T-I- (x,7) 

= T- (x,a) T- (x, 7 ) |r+ {x,a),T + (x, 7 ) )■ -f T_ (x,Q;),r_ (x, 7 ) \f+{x,a)f+ (x,j) 

+ 


1 2 

r(a,j),T (x,a) T (x, 7 ) 


12 2 
r (x, a), T- (x, j)r(a,-j + S) T+ (x, 7 ) 


-f 


2 1 1 
r (x, 7 ), T- (x, a)r(-a + S,j)T+ (x, a) 
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After taking trace on T, we find 


1 2 
tr T (x, a), tr T (x, 7) ? = tritr 2 


T (x,a),r (x, 7 ) } = 0, 


that is 

{r (a),T (7)| = 0. (19) 

In other words, T(a) constructs a one-parameter involutive family. Here we remark that no symmetry 
conditions of r- matrix is used to obtain equation dil)- Consequently it can be applied to general model. 


3.2 Conservation condition 


If T{a) is a generating function for integrals of motion, it must be independent of time. We find 
dt trT (x, a) 

= dt tr a)T(x+,x-,t, a)K-{t, a)T“^(x+, x_, t, —a -I- ^)} 

= tr {[dtK^{t, a) — V (x+, t, —a + d)K^(t, a) + K^(t, a)V (x+, t, a)] 
xT(x+, x_, t, a)K-{t, a)T~^{x+,X-,t, —a + S) 

+ [dtK-(t, a) — V(x-,t, a)K-(t, a) + K-(t, a)V (x_, t, —a + J)] 

xT“^(x+, x_, t, —a + 6)K^{t, a)T(x+, x_, t, a)} . (20) 

Taking dttr'T (x, a) = 0, and supposing there has no connection between boundary variances on each end, 
we obtain the evolution equations of K± matrices 

dtK+{t,a) — V{x+,t,—a + S)K+{t,a) + K+{t,a)V{x+,t,a) = 0, 

dtK-{t,a) — V{x-,t^a)K-(t,a) + K-{t,a)V{x-,t,—a + S) = 0. (21) 


For these equations, we find immediately that there have two isomorphisms between and AT-, which 
are Ar+(a) ^ K-(—a + S) and iF+(a) ^ KZ^{a:). 

If K± are constant matrices {dtK± = 0), the equation ( pl|) can be simplified as 

V{x+,t,—a + S)K+{t,a) = K+(t,a)V{x+,t,a), 

V{x-,t,a)K-{t,a) = K-{t, a)V{x-,t,—a + 6). (22) 


In the present case, we remark that K± matrices are not singular matrices, so that determinants of V satisfy 


detV{x±,t, —a -I- ^) = detV{x±,t, a), 


(23) 


by which we can obtain the value of 6. After inserting <5 value into equation (22), we can find some nontrival 


CIBC when nontrivial K± matrices appear. In other words, a class of K± matrices is relate to a class of 
integrable boundary condition. 

We must point out that K± matrices depending on field variables take its meaning in fact [§. On one 
hand, we must use such K± matrices in order that quantity (0) can be regarded as a generating function on 
periodic boundary condition too. On the other hand, studying such K± matrices, we understand integrable 
condition more deeply. 

When Sklyanin’s function (^ is regarded as a generating function in ATFT, boundary K± matrices will 
have no constant solution in equation dil) (except for Sine-Gordon theory). So we have to solve equation 
( pT|) as differential equations. Besides these difficulty, even when one had found a nontrivial solution, he 
should have to prove the involution condition again because condition (0 may be broken. In our method, 
6 added on spectral parameter guarantees the existence of constant K± matrices, and d can be solved by 
means of equation (|^), so K± matrices solving procedure is simplified effectively. 


Proposition 3. If K± matrices satisfy not only algebra equation (0) and (^^, but also evolution 
equations of trT (x, a) is a generating function for integrals of motion. 
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There has a difficult step in proving this proposition. Are solutions of K± matrices in equation 
compatible with its algebra equation (|^) and ( 0 ) ? Although one found it is true in Sine-Gordon theory 
and it has been proved in ATFT^j with the form of ([T^), it still keeps an open problem in general 
theory. If this compatibility is satisfied, the proposition is proved naturally. 

3.3 Other generating functions 

As exhibited in above subsections, we also obtain other K± matrices’ algebra and evolution equations when 
quantity (@) or is regarded as generating function for integrals of motion. In the form of (p^, it read 
as 

0 = + 6,-^ + 5)K+{t,a)K+{t,"i) + + 5)K+{t,"i) 

+K+{t, {-a + (5, -i)K+{t, a) + K+{t, a)K+{t, 7 )r(a, 7 ), (24) 

0 = r{a,j) K- {t,a) K- K- +6,"f) K- 

+ K- {t, 7 )r‘= (a, -7 + S) K- {t, q;)+ K- (t, a) K- {t, 7 )r*i‘^ {-a + <5, -7 + <5). 

where the upper index ”ti, i = I, 2” denote transposition on the space, its evolution equations read as 

0 = dtK+{t,a) + V*{x+,t,—a + S)K+{t,a) + K+{t,a)V{x+,t,a), 

0 = dtK-{t,a) — V{x-,t, a)K-{t, a) — K-(t,a)V*{x-,t,—a + 6). (25) 

Let dtK± = 0, we can also obtain 6 value in ( ^^ by taking determinants. With constraint conditions of 
( p4| ) and (p^), quantity (^ is a generating function for integrals of motion. 

In the form of (]l^, the similar equations are still balance except that we must use ’’f” (conjugation) 
instead of (transposition). If those modified equations are satisfied, quantity 0 ) can be regarded as a 
generating function for integrals of motion too. 

Now, we have obtained three forms of generating function as well as their constraint conditions. If all 
of them are regarded as generating functions, we believe that they are the same one in fact. In the next 
section, we will prove it explicitly in ATFT. 


4 Classically Integrable boundary condition in ATFT 

4.1 Links among generating functions 

Lagrangian in ATFT on independent boundary condition is §0: 


+ 00 +00 

C = J dx J dt 

— 00 —00 


{ 1 2 

0{x - x-)0{x+ - x)[-d^(j)adf^(l3a “ ^ ^ ni(e^“‘ ‘^ - 1)] 

-6{x - x-)V-{(j){x-), a^0(x_)) - 5{x+ - x)V+{(j){x+), d^,(j){x+))} , 


(26) 


where m is mass scale and P is the coupling constant in real domain; at are simple roots of a simple Lie 

r 

algebra of rank r (included the affine root ao). There are '^mai = 0 and no = 1. It is a theory of r scalar 

0 

r —1 

fields {ai ■ (j) = The potentials V+ and V- are additions on the ends x+ and x-, respectively. 

a—0 

They denote independent boundary conditions. When V± depend on (l){x±) only (and independent of its 
derivatives), we obtain: 


{d1-dl)(j) = x-<x<x+, 
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dV± 


X = x±. 


(27) 




Lax pair in ATFT read as (A = e“) 

Uix,t,X) = - |i/3i7 + 

V{x,t,\) = - |i/3i7 • i93,0 + m^7r^(A£;ai - . (28) 

in which H and E±on are the Cartan subalgebra and the generators responding to the simple roots, respec¬ 
tively, of the simple Lie algebra of rank r. The coefficients rm are equal to ruaflS. There are the Lie algebra 
relation: 

[Hi,Hj\={) , [H,E±a^] = ±aiE±a^, 

[E^„E_^,]=2a,-Hl{a}). (29) 

It is pointed out by Hollowood 111 that the complex affine Toda theories have soliton solution (in 
which coupling constant /3 is purely imaginary), in contrast with the real coupling constant ones. And its 
Lagrangian, motion equation and Lax pair can be expressed similarly to equation (^)-(^) except for taking 
(3 ^ i/3 G Re). In our paper, we use (^)-(^) equations in general, and distinguish them only when the 
real and imaginary cases can’t be treated in the same way. 

For those generators in Lax pair (p^, we can find a representation in which they satisfy: 

Hl = H\ = H, ,Ei^^=Ei^^=E^^,. 

So there has an automorphism map : 

H, H'i = = -Hi , 

Eon *■ OLi = Ect^^ = E-on , 

E—cti —*■ E' — a- = E—ctf^^ = Eon ■ 

The new generators satisfy the same Lie algebra relation (|2^). In other words, there are 

U\x,\) = U{x,\-^) = -Vl-^U{x,-X)n , 

V\x, X) = V{x, -A-i) = -Xl-^V{x,-X)n . (30) 

We remark that equation ( |30| ) can be applied to both real and imaginary coupling constant cases. But if 
one uses ’’f” instead of ”f’, equation (^ ) must be modified because of its complex fields. 

From definition of T{x, y, t, A), there has 

dxT\x,y,X) = T\x,y,X)U\x,X) 

= T\x,y,X)[-n-^U{x-X)Vi\, 

or 

d, [L!T*(cr, 2 /,A)II-i] = - [lIT*(x,y. A)!!"!] U{x,-X). 

Comparing it with dxT~^{x,y, X) = —T~^{x,y,X)U{x,X) and the initial condition in (||), we obtain 

nT\x, y, X)n~^ = T"^(x, y, -A), 
or 

VtT^{x, y, a)n“^ = T~^{x, y,a + in). (31) 

It means 


tr{K-(a)T ^{—aX-5)Ko^{a)T{a)'\ 

= tr {K-{a)nT*^{—a -I- <5 -I- f7r)n“^Ar+(a)T(Q!)} 

= tr^K_ia)T\-a + S')K+{a)T{a)Y (32) 


7 



in which K± matrices in ( 0 ) and ( 0 ) are distinguished by K± and K± now, and the quantities added on 
spectral parameter become <5 and 6' respectively. In other words, quantity & is equal to (0), if 6' is equal 
to ^ + *7r and reflection matrices satisfy 

K_{a) = K_ia)n , K+{a) = n-^K+{a). (33) 


Using the second equation of ( p0| ) and comparing equation (0 with we find these relations appear 
again. So quantities (0 and ( [l^ are the same one in fact, when both of them are regarded as generating 
functions for integrals of motion. 

In real coupling constant and real fields case, if we use ’’f” instead of ”t”, equation ( ^0|) is still balance 
when spectral parameter is real. So we obtain a relation similar to again. In this case, when we rewrite 
as T(a) = trK-{a)T^—a + 6'')k+{a)T{a), then 


K±{a) = K±{a) , 6' = 6" = S + iir. 


(34) 


Now, we have proved quantity (0 is equal to 0 ) when both of them are regarded as generating 
functions in ATFT. When coupling constant is real, they are equal to generating function (|^ too. But 
when coupling constant is purely imaginary, equation ( ^0[) may be not satisfied, so K± matrices in (^) may 
have no constant solution. 


4.2 Classically Integrable Boundary Condition 


In real coupling constant ATFT, if we regard T(a) = trK-{a)T^—a+6)K+{a)T{a) as a generating function, 
we will obtain the evolution equation of K+ in a;+ boundary: 

U^(a:+, —a + 6)K+(a) + K+{a)V (a;+, a) = 0. (35) 

After taking 5 = 0, we obtain equation of .A^^(A)(in which A = e“): 


1 

2 



m 


-I + 


i?+i(A), ^ 

0 


(36) 


By boundary equation (p^, it is just the reflection equation appears in paper Q. We find K+ and T 
matrices defined in that paper are just the quantities of and T~^ in our paper, according to different 
definition of Lax pair. Analogy with method in paper we solve equation (^) and obtain CIBC in ATFT. 
In simple-laced case, it is the same as paper [|| 


—0x9 

m 


in which \Bi\ 
or 


r 


-Y.B, 





2 = 0, l,---,r, 

0 = 0 , 


(37) 


In imaginary coupling constant one, we regard T{a) = trK-{a)T*{—a + S)K+{a)T{a) as a generating 
function. The results in real coupling constant can be used, thanks to section 4.1. In other words, the new 
CIBC can be obtained by analytic continuation by /3 ^ f/3 (/3 G Re). It reads: 



m 

in which \Bi\ 
or Bi 


r 







2 = 0, l,---,r, 

0 ,J = 0, l,---,r. 


(38) 


We remark that Sklyanin’s method ||^ can’t be used in ATFT except for Sine-Gorden theory, this con¬ 
clusion comes from the fact that 5 0 on independent boundary condition if we regard 0 ) as a generating 

function. Now, we must take S = —in according to equation ( ^ ) and (p^. Sine-Cordon theory is an 
exception in which it is satisfied both 5 = 0 and 5 = —in. 

As we discussed in section 3, in the classical case, no symmetry conditions of r- matrix is necessary in 
constructing generating function for integrals of motion. So it is interesting to study whether commutative 
quantities can be constructed with less symmetry of R- matrix in the quantum case. 












5 Quantum integrable systems with boundary 


There are many papers (for example, and [p^ - fl^ ) in which the authors deal with quantum integrable 
boundary condition in two-dimensional lattice models. As far as we know, both unitarity and crossing 
unitarity conditions (or the weaker property fl^ ) are used in constructing commutative quantities. Since 
finding crossing unitarity condition of a given R- matrix is a difficult problem, it is useful to construct 
commutative quantities without this symmetry. 

In this section, we explore how to obtain commutative quantities by means of unitarity condition only. 
Unitarity condition read as: 

i?i2(n)i?2i(-u) = ^(u), (39) 

where ^(u) is some even scalar function and R- matrix is a solution of quantum Yang-Baxter equation 
(YBE): 

Ri 2 {u - v)Ri3{u)R23{v) = R23{v)Ri3{u)Ri2{u - v). 

As usual, the transfer matrix t(u) is defined as 

t{u) = tr Tj^{u)T-{u), (40) 

and each entry of 7i|_(M) commutes with T-{u). 


Proposition 4. If T± satisfy such equations 


in 




in 


T + {u)R*^^{u+ -5)r+ {v) = T+ (?;)i?i" 2 (^^-H - ^) T + (n)i?i 2 (-M-), 
Ri2{U-) T- (u)i?i2(“^-l- + (^) = (^)^21*^(^-) 


(41) 


and the quantum R- matrix obeys the unitarity condition, then, transfer matrix t{u) defines a one-parameter 
commutative family. 


For explicity, we use ^i ^ and ^2 ^ replace ^ ^(■u+ — (5) and ^ ^(—it_) respectively, as well a.s u± = u±v. 
The proof is directly: 

t{u)t{v) = triT+ (u) T- (u) tr 2 T + {v) T- {v) 


tri 2 T+ (m) T+ {v) r_ (u) T- (v) 


in 2 


I ‘1 


= tri2 T+T+ R2 i{u+ - 5)Ri2[-u+ + S) T_T- 

= tri2{T+ Rl\iu+ - 6) tIYHT- R[\{-u+ + 6) 

= tri2{T+ R2 \{u+ - S) T+Y^*^R2i{-u-)Ri2{u-){T- RlY-u+ + 5) T-} 

= Cr'? 2 “' tri 2 {i?*f (-n_) r+ Rl\iu+ - S) tIY^^HRi 2 {u-) T- R[Y-u+ + S) T-}, 
using equations (|4l|), we find 


A 2 


in 


tiu)t{v) = ^ 1 - 1 ^ 2 “' tri 2 {T+ R[Yu+ - S) r+ Ri 2 i-u-)Y^*HT- Rl\{-u+ + S) T- 


= 4" tn2{T+ RlYu+ - S) T+Y^*HT- Rl\{-u+ + <5) T-} 




= tri2{T+ RlYu+ - 5 ) T+YYT- R 2 Y-U+ + S) T-} 


+ -“-12 
2*2 1 2*2 1 

tri2 T+T+T_T- 


*2 
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In quantum spin chain model, it is convenient that T± take such representations 


T+{u) = K+{u), 

T_{u) = T{u)K_{u)T\-u + 5) (42) 

= L]\f{u) ■ ■ ■ L2 {u)Li{u)K_{u)L\{—u + 5)L\{—u + ^) • • • + 6), 

in which transposition acts on the auxiliary space and n = 1, 2, N denote quantum space. There has 
a relation between R and L operators 

Rabiu - v)La{u)Lb{v) = Lb{v)La{u)Rab{u - v). (43) 

Let T-{u) = Lpf{u)T_!_{u)L*jy{—u + d) and insert equation ( ^ ) into equation (^l|), we find the second 
equation of (|l]) becomes 

l.h.s. = Rab{u-)LaN{u) RL + <5)i?‘“f,(-u+ + 6)LbN{v) TL + S) 

= Rab(U-)LaN(u) T- LbN{v)R^^bi-U+ + 

= LbN{v)LaN{u)Rab{u-) T!_ i?*'),(-M+ + 5) Ti + S)LI’’j^{-U + (5), 


r.h.s. = LbN{v) TL + 5 )RI^^{-u+ + 5 )LaN{u) T( L^“^(-u + ( 5 )i?j“‘^(u_) 

= LbN{v) r LMu)Rl'’a{-u+ + 5 )LI\{-u + 5 ) r L^“^(-u + S)Rll^>’{u.) 
= LbN{v)LaN(u) TL Rba{-'^+ + <^) + ^)LbNi-U + ^)- 


In other words, this equation is reduced to 


Rabiu.) r Rlii-u+ + 5) r=r Rti-u+ + s) r Rlt^iu.). 

We proceed to do the above reduction repeatedly until all of L operators beside K. matrix disappear. 
At last, we obtain the reflection equation about K. only. Now, reflection equations of K± are 


I *2 


. t2 


R-ii'^i-U-) K+ {u)R^^iu+ - S) K+ (v) = K+ iv)RL2iu+ - S) K+ {u)Ri2{-U-), 

Ri2{u-) K- (u)i?j2(“^+ + ^) K- i^) = K- i'^)R*2ii~'^+ K- i'^)R*2i^i'^-)^ 

and transfer matrix t(u) becomes 

tiu) = tr Ar+(u)T(tt)Ar_(u)T*(—u + <5) . 


(44) 


(45) 


We remark that there has no obvious relation between K± matrices. If some symmetry conditions are 
used, relation between and K. matrices may be found. 

For example, R- matrix in paper [O has PT symmetry and crossing unitarity 


i?12(u) = i?‘\‘^(u), 

1 , 1“^ 

Ri 2 {u) = VRili-u-p)V ■ 

By PT symmetry, we find there has an isomorphism between boundary matrices: 

K-iu) = K^j^i—u + (5). 

If both PT symmetry and crossing unitarity are considered, there is another relation as 

K_iu) = K-^iu +p)M-^, M = V*V. 


(46) 


(47) 


(48) 
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From and (|4^), it means 


(49) 


K*^_(—u + 6)MK^{u + p) = Id, K*_{—u + S)K_{u — p)M = Id. 

These equations may be regarded as constraint conditions on 5. 

Now, it is interesting to compare our commutative quantities with those of Mezincescu and Nepomechie 
[ p^ . Using conditions of (^) and unitarity condition i?i 2 (u)i? 2 i(—u) = ^(u), we obtain 

R[],i-u-p) = ivY-R%^^{u){v 'y- 

= ^{-u){vY^ R-Yi-u)iv 'y\ 


or 


Ri 2 (—u + i5) — —u p 6) (U)*^ -^ 12 ^ (—^ P d) {V )*^ • 

If Ln{u) is defined as Ln{u) = Lan{u) = Ran{u), we obtain 

= ^{—u + p + 5)V*'L^^{—u + p + 5){y 
T*(—m + 5) = L\{—u + 5)L\{—u + 5)---V'j^{—u + 5) 

= i^[-u + p + 6)V*T-\-u + p + 5){y-^Y. 


In other words, transfer matrix 


becomes 


t{u) = tr K+{u)T{u)K_{u) +p + 5)V^T-^{-u +p + 5){V-Y^) 

= I,^{-u +p + 5)tr {{V-^YK+{u))T{u){K_{u)VYT-^{-u +p + 5). (50) 

It is convenient to multiply (|50| ) by {—u + p + (5) before we regard it as the commutative quantities, 
i. e., 

t{u) = tr KY{u)T{u)K'_{u)T~^{—u + p + d). 

If d is equal to —p, it is just the one in paper 

As one of the main results in our paper, we have constructed the commutative quantities with unitarity 
condition of quantum R- matrix only. As discussed in the classical case, with symmetry conditions of (^^, 
we can find another form of commutative quantities, and these two forms are the same one in fact when 
both of them are regarded as commutative quantities. 

Finally, we study the classical counterparts of reflection equation s (4^ ) by modifying unitarity condition 
to Ri 2 {u)R 2 i{—u) = Id. In the classical limit, as ?i ^ 0, one has 


[ , ] = —ih{ , } ; R{u) = Id + ih r(u) + o[TY). 


So unitarity condition means ri 2 (u) = —r 2 i (—u) and quantum YBE goes over into the classical YBE. We find 
reflection equations ( p4[ ) just turn into equation ( p4[ ) in which r{a,P) is equal to r{a — (3) now. In contrast 
with the quantum case, there has an isomorphism between and K-, which is A'_|_(a) —> KZ^{cx). 


6 Conclusion and discussion 


In this paper, we obtain three possible generating functions for integrals of motion in classically integrable 
field theory on a finite interval with independent boundary condition on each end. As constraint conditions, 
we find K± matrices’ algebra and evolution equations. In contrast with other’s methods, a new parameter is 
added on spectral parameter, and we expect it shall simplify the procedure of solving K± matrices effectively. 
In ATFT, we prove these generating functions are equivalent to each other and its links are discussed too. 
Our results show that two of these generating functions are always valid in both real and imaginary coupling 
constant cases. 

It is remarkable that no symmetry condition of r- matrix is used when we regard quantities (0), ® 
and 0) as generating functions for integrals of motion, so we expect it shall be applied to more integrable 


II 




models than I, I- As demonstrated in section 4, the added parameter 6 improves this possibility. 

We also extend our results to quantum spin chain, we have proved that unitarity condition of quantum 
R- matrix is sufficient to construct commutative quantities with boundary. Reflect equations of K± are 
obtained. Relation between boundary K± matrices is found when PT symmetry and crossing unitarity 
condition of R- matrix are considered. With these symmetry conditions, we also find another form of 
commutative quantities as the one defined in [p^ . Finally, we find that classical counterparts of quantum 
reflection equations is just the one which we obtain by classical quantity 
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